Neutrino spectrum from the pair-annihilation process in the hot stellar plasma 
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An new method of calculating the energy spectrum of neutrinos and antineutrinos produced in 
the electron-positron annihilation processes in hot stellar plasma is presented. Detection of these 
neutrinos, produced copiously in the presupernova which is evolutionary advanced neutrino-cooled 
star, may serve in future as a trigger of pre-collapse early warning system. Also, observation of 
neutrinos will probe final stages of thermonuclear burning in the presupernova. 

The spectra obtained with the new method are compared to Monte Carlo simulations. To achieve 
high accuracy in the energy range of interest, determined by neutrino detector thresholds, differential 
cross-section for production of the antineutrino, previously unknown in an explicit form, is calculated 
as a function of energy in the plasma rest frame. Neutrino spectrum is obtained as a 3-dimensional 
integral, computed with the use of the Cuhre algorithm of at least 5% accuracy. Formulae for the 
mean neutrino energy and its dispersion are given as a combination of Fermi-Dirac integrals. Also, 
useful analytical approximations of the whole spectrum are shown. 
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I. INTRODUCTION & MOTIVATION 

We present here a new method which allows us to cal- 
culate the spectrum of neutrinos produced by thermal 
pair-annihilation processes in the hot plasma in the core 
of massive pre-supernova star. The paper is organized 
as follows: In sect. II the electron-positron pair annihila- 
tion processes into neutrinos in the Standard Model are 
briefly summarized. Monte Carlo simulations are dis- 
cussed in section All Al Simple estimates of the average 
energy of neutrinos based on the Dicus cross-section [jj 
are presented in Subsection IIII Bl 

The complete spectrum is obtained using the differen- 
tial cross-section (Eq. derived in Subsection IIII CI 
where the spectrum is given in the form of three dimen- 
sional integral, that is easy to evaluate numerically with 
the Cuhre or Monte Carlo algorithms. Additionally, we 
express the moments of the spectrum as a combination of 
Fermi-Dirac integrals (Subsection llll D() . that are used to 
obtain convenient analytical approximation of the spec- 
trum in Subsect. llIIEl 

Our long-term goal is to explore the entire neutrino 
spectrum produced by massive stars at late nuclear burn- 
ing stages. Unlike the Sun, these stars after carbon ig- 
nition emit both neutrinos from weak nuclear reactions 
and thermal neutrinos. Moreover, photon luminosity is 
negligible compared to neutrino luminosity and therefore 
these objects are referred to as neutrino-cooled stars 0- 
Neutrinos emitted as neutrino-antineutrino pairs domi- 
nate up to the silicon burning |3( ■ The electron antineu- 
trinos are much easier to detect than neutrinos 0,|HSI3- 
We consider them first. 



Production of pair-annihilation neutrinos is the dom- 
inant thermal process with relatively high average neu- 
trino energy of (£„) ~1.5 MeV. However, the informa- 
tion on neutrino emission is not complete without full 
knowledge of all neutrino processes. This is important 
because other neutrino processes can also produce low 
energy neutrinos which are difficult to detect. These neu- 
trinos somehow "steal" detectable energy from the stellar 
core reducing chances for e.g. supernova prediction. 

The results presented here are general in nature, 
and are valid for pair-annihilation neutrinos from e + e~ 
plasma in the full range of temperature and chemical po- 
tential. 



II. PAIR ANNIHILATION IN THE STANDARD 
MODEL 



According to the Standard Model of electroweak in- 
teractions, the electron-positron pair may annihilate not 
only into photons but also into neutrino-antineutrino 
pair: 



(1) 



In the first order calculations, sufficient in the consid- 
ered energy range of several MeV (E <C M w ± z <> where 
^w ± ,z° is the intermediate boson mass of ^100 GeV), 
two Feynmann diagrams (Fig. ^) contribute to the anni- 
hilation amplitude. 

The contribution from the left diagram in Fig. ^ i.e. 
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FIG. 1: Feynmann diagrams leading to the first-order ampli- 
tude for e + e~ annihilation into neutrinos. Electron neutrinos 
are produced via both diagrams, while p and r neutrinos only 
by " Z° decay". 



boson exchange (charged current) is given by Q : 
iGp 



V2 



U„j a (l - 75)u e V e ^a(l ~ JB)V V . (2) 



This part of the process produces only electron neutrinos. 
Using the Fierz transformation: 

a[7 M (l-75)]6 c[ 7/1 (l- 75 )]d = -o[y*(l-7fc)]d c[ 7m (1- 75 )]6 

(3) 

we may write @ as 1 : 
iG F 



V2 



W„ 7 "(l - 7 5 )v„ We7a(l - 75)Me 



(4) 



The Z° boson exchange (neutral current) gives: 
iG F 



V2 



u^ a (l - 7s )i>„ v e j a (g v - g A To)u e (5) 



where: 



1 , 1 

9v = + 2sm w, 9A = -^ 

By adding jj to O we obtain the total annihilation 
amplitude: 

_A4 — ^ F -. _ r, r« err/ _ r</ 
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Wiy7 Q (l - 75)«i/ V e ^a{Cy - C' A J5)u e - (6) 



Here Gj? = 1.16637(1) x 1CT 5 GeV~ 2 is the Fermi con- 
stant of weak interactions [t| . The parameters C v and 
C A (/ = e, fx,T - neutrino flavor) in the Standard Model 
of the electroweak interactions are: 
for electron neutrinos: 



C e v = I + 2 sin 2 



(7) 



for \i and r neutrinos: 
1 



C<p T = -~ + 2sm 2 e w , W4 



(8) 



where 9w, the Weinberg angle, is sin 2 6w — 0.22280 ± 
0.00035 0. 

As the pair-annihilation into electron neutrinos pro- 
ceeds via both charged and neutral currents we have: 



CI 



C?/ T + 1, 



c 



c>i' T + i 



A-^A I"- (9) 



Annihilation into fi and r neutrinos proceeds via the neu- 
tral current only. 

In general, the entire spin-averaged information 2 on 
e + e~ annihilation process Q is included in the spin- 
averaged (<7i - spin states) squared annihilation matrix 
M 2 : 



(10) 



(7 J- ,(T C7 V ,(Jn 



The amplitude M 2 can be calculated using e.g. the 
Casimir trick: 



M 2 = 8G F 2 { U 

(< 



C{ - C; 



Pi ■ Qi Pi ■ Q 2 



+ (c f A + Cfj P 2 ■ Qi Pi ■ Q 2 

+ m e 2 {c v 2 - C f A 2 ) Qi-Qa} (11) 

where Pi = (-Ei,pi), P 2 = {E 2 ,V 2 ), Qi = (£i,qi) 
Q2 = (£2^2) is, respectively, the four- momentum of the 
electron, positron, neutrino and antineutrino. 

According to the definition of the cross-section da, the 
number of collisions dN occurring in volume dV in time 
<ftis[n|: 



dN 
dVdt 



= dav dni dm, 



(12) 



where dni y2 are particle densities. In the case of two 
incoming particles, the differential cross-section can be 
computed from the Fermi Golden Rule formula and is 
given by the expression : 



d<TV ~ 2E X 2E 2 (2tt) 2 *' 



6*(Pi+P 2 -Qi-Q 2 ) M 2 



2Si 2£ 2 



(13) 

Hence, if M 2 is known, the total neutrino emissivity 
from e + e~ plasma can be calculated by performing ap- 
propriate integrations |l3[ : 



rf 3 pi d 3 p 2 rf 3 qi d 3 q 2 
2E\ 2E 2 2E\ 2£ 2 



AA/2<5 4 (P! + ...) M 2 



(14) 



1 We are dealing here with the field operators, not numerical 
spinors, so the overall minus sign in Eq. disappears 0. 



2 If electrons are polarized then we have to calculate Af 2 using the 
density matrix for polarized fermions llOl . 
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Substituting into the integrand of i|14|) the appropriate 
expression for the factor A we obtain formulae for the 
number emissivity (A = 2), the total emissivity (the 
energy carried by neutrinos and antincutrinos) (A = 
E1+E2 =81 + 82), the antineutrino emissivity (A = £%) 
etc. Functions /i,/2 are the Fermi-Dirac distributions 
for electrons and positrons, respectively: 

ft = gCBi-^/feT + P h = e {Ez+n)/kT + I' ( 15 ) 

where /i is the electron chemical potential (including the 
rest mass) and neutrinos and antineutrinos are assumed 



to escape freely from the plasma. 

The total neutrino emissivity, as required e.g. for stel- 
lar evolution codes, including the emission of all three 
flavors (e, /x, r) is the sum of three terms of the form of 
the integral (|14(l with appropriate coefficients (0 IBJl in 
CD- 

The formal expression 114|l must be transformed into 
analytical or convenient numerical form before actual ap- 
plication for e.g. stellar evolution codes or signal estima- 
tions in neutrino detectors. The integral 114J above is 
simplified significantly if Lenard's formula |l4l Il5|: 



IsT^sT Q?Q * 54(Pl + P2 - Ql - ° 2) = Ya [ aaP{Pl + P2)2 + 2 (jP " + p ^ + p ^} 6 + P2)2 ] (16) 



is applied to integrate over outcoming neutrino momenta giving the well-known cross section 3 of Dicus for 
annihilation of the e~e + pair with four-momenta Pi and P2: 



<7V D 



12tt ExE 2 



1 + 2 



Pi -Pi 



(17) 



The formula l|14|l can be now transformed using i|17|) 
into: 

J d 3 pid 3 p 2 A av D /i/a. (18) 

In contrast to a more general Eq. 1|14(1 . in (|18|1 A must 
be a function of e~, e + momenta pi, P2 only, and we are 
unable to compute the antineutrino spectrum separately. 
The expression (|18|) above is actually a three-dimensional 
integral - due to the rotational symmetry only lengths 
and angle between pi and P2 are independent variables 
of integration. Therefore, (|18|) may be integrated, leading 
to a combination of Fermi-Dirac integrals (cf. Sect. lIII B|l 
which are easily computed numerically. 

We can obtain the total neutrino luminosity in this 
way. Unfortunately, if one attempts to consider the de- 
tection of these neutrinos very important information 
about the neutrino energy is missing. Detection efficiency 
depends strongly on the neutrino energy - higher energy 
neutrinos are much easier to detect - so the spectrum 
parameters and the shape actually decide whether these 
neutrinos will be observed or not. Therefore we must go 
back to the general formula (|14f) and try to attack it in 
an another way. We consider various possible approaches 
in the next subsections. 

3 We denote expression 1171 as crvjj to avoid confusion with dif- 
ferential cross-sections 1311 . 



III. NEUTRINO SPECTRUM 
A. Monte Carlo simulations 

The Monte Carlo simulation is a method often used to 
calculate annihilation spectra. We applied it in our pre- 
vious work 0, with some minor errors in the algorithm. 
To improve the calculations we develop semi-analytical 
methods, to make cross-check of the obtained results. 
Poor accuracy for higher energies of previously obtained 
results did not allow us to compute reliably the spectrum 
for energies E>4 MeV. Using better random number gen- 
erator, namely Mersenne Twister routine [16|, and longer 
simulation runs, we were able to extend the spectrum to 
6-7 MeV at least, but accuracy was still unacceptable, 
cf. Fig.EJ 

As we can see in Fig. the Monte Carlo simulation re- 
produces the spectrum well. Unfortunately, this method 
is very slowly convergent: time required to compute the 
grid of neutrino spectra for e.g. the purpose of inte- 
grating them over the entire volume of the hot stellar 
core is very long. Moreover, the accuracy of this method 
is limited by resolution of the random number genera- 
tor, and errors in the interesting range above ~5 MeV 
(Super-Kamiokande threshold) are large. Therefore de- 
velopment of the analytical formulae is welcome. Never- 
theless, Monte Carlo simulations based only on the am- 
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FIG. 2: Spectrum resulting from the Monte Carlo simulation 
for 10 6 , 10 7 , 10 s , 10 9 , 10 10 and 10 11 loop size (symbols) versus 
spectrum from 137H computed with aid of Cuhre algorithm 
(solid line, relative accuracy 0.05). We can notice decreas- 
ing errors of the Monte Carlo simulation near the maximum 
and increasing number of events in the high energy tail. Obvi- 
ously, the Monte Carlo simulation is unable to reach accuracy 
of Cuhre methods for E 3>5 MeV. This figure clearly justifies 
efforts to compute the spectrum using the appropriate cross- 
section 1311 1. instead of using Monte Carlo simulation based on 
the knowledge of the annihilation matrix 11 II . We would like 
to point out that three-dimensional integral 1371 may also be 
computed using advanced Monte Carlo adaptive algorithms, 
giving similar spectra much faster, but they may fail in some 
cases, cf. Fig. [3] Spectrum computed for typical situation 
in stellar center during Si burning with T = 0.319 MeV and 
H = 0.85 + m e MeV, 



plitude M 2 flll|) may be used even in the case where 
analytical formula for the cross-section is not available. 



B. Emissivity and average neutrino energy 

Some quantities related to the neutrino emission pro- 
cess, such as reaction rates, the total emissivity and the 
average neutrino energy, can be easily obtained directly 
from the Dicus cross-section (|17[) as electron energy mo- 
ments. These quantities are average values for v — v 
pair. However, in Sect. IIII CI we show how to modify 
average values by small additive terms to obtain the 
relevant quantities for neutrinos and antineutrinos sep- 
arately. This becomes possible only using a more general 
cross-section <|31[1 . 

Our intermediate goal is to compute the integral: 



In 



(2n) 



d d Pl d 3 p 2 (jv e) E 1 n E 2 m fih 



(19) 



i.e. the expression {TBI) with A = Ei n E 2 m where the 
cross-section avr> is given by the formula (|17|1 . and n,m 
are arbitrary real numbers. 



The integral over cosine of the angle between pi and 
P2 (cos 9) can be easily done, while integrals over dE^dE 2 
can be separated and expressed by the products of Fermi 
integrals. This is a result of the fortunate coincidence: 
the expressions \/E 2 — m e 2 appear (i) from phase-space 
factors d 3 p and (ii) from four-momentum products, but 
the latter ones are always squared because the integration 
of cos™ 9 with an odd power gives zero. 

Let us define functions 



M nm _ (j Cv 2 2Cj 2^ GT /2 _ i/2 G± /2 _ i/2 



9Cy 2 G;f /9 G; 



n/2^ m /2 + ( C V H 
^71/2-1/2" m/2+1/2 



Ca) ( 4G ^/2+l/2 G m/2+l/2 
" < ^n/2+l/2 < ^m/2-l/2) ( 2 ^) 

where the Fermi integrals G^ are: 
1 



v 3+2n 



1 + cxp(x ± (3) 



dx 



(21) 



a = m e /kT, (3 = jjL e /kT, x — E/kT. 
The functions M™ 1 obey the relation: 

M nm _ M rnn 



(22) 



For example, n-th moment of the electron energy is 
proportional to M™° and n-th moment of the positron 
energy is proportional to M™° = M°™ 

Electron (positron) energy moments can thus be cal- 
culated as: 



8+n+m 



18tt 5 



-M r ' 



(23) 



We can now write elegant expressions for the emissiv- 
ities. which are equivalent to well-known formulae |17| . 
0, 0|- The total emissivity (the neutrino energy pro- 
duced per unit volume and unit time) is equal to: 



dE 



M 



10 \ 



(24) 



dVdt ~° 18tt 5 
while the number emissivity (particle production rate) is: 



dN 
dVdt 



F = 2R 



18tt 5 



(25) 



The reaction rate R is half of the particle emissivity as 
two neutrinos are produced in a single event Q). In 
eq. H24|) energies of two neutrinos [y and D) are added 
together. 

We can also compute the average neutrino energy from 
the pair-annihilation process: 



Q = Q_ 

F 2R 



(26) 
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Average neutrino energy as a fraction of the electron 
rest-energy m e c 2 is: 



Non-relativistic and non-degenerate plasma 
kT <C m e , kT > [i: 



(£) = 



M 



Ml 



M oo + M oo 



\m e c 



(27) 



(£) = m e + - feT 



(29) 



Unfortunately, no more informations of interest on the 
neutrino spectrum can be extracted from Lenard-formula 
based calculations. We wish to point out that the equa- 
tion (|2*7|) gives the average v x -v x energy (£) = {{£\j + 
(£ 2 })/2. However, both mean energies (£\) and (£2) are 
not identical, but differ slightly. They can be both com- 
puted from eq. i|39|) . 

Knowledge of the mean neutrino energy from the for- 
mula (|27(l allows one to quickly estimate chances for neu- 
trino detection from a given reaction. At present, neutri- 
nos of a few keV energy are impossible to detect, 1 MeV 
neutrinos are difficult, while for 10 MeV neutrinos we 
have mature detection techniques with megaton targets 
soon available. 

Using analytical approximations for Fermi-Dirac inte- 
grals [jjj, |2(j, we derive useful analytical formulae for 
average neutrino energy. Particularly simple expressions 
exist in the following cases: 

• Relativistic and non-degenerate plasma kT > 2 m e , 
kT > fi: 

= 2700 C(5) j fcT = 4 10 6fcT~4jfcT (28) 

7 7T 4 



• Degenerate case fi 3> kT: 



, , 2 m e 
(£) = -^ + 2kT+^- 



(30) 



Unfortunately, in the central region of a massive star 
m e ~ kT ~ fi and none of the cases above holds. Eq. 127|) 
must be evaluated numerically. 



C. Spectrum shape 

To find the flux of more-than-average energy neutrinos 
we have to compute the neutrino spectrum. To do this we 
go back to the formula l|14|) , and integrate it step-by-step. 
Without the Lenard's formula H16fl these calculations are 
"tedious algebra" [2l| , j2^| and lead to very complicated 
expressions. 

We have computed the essential differential cross sec- 
tion for e + e~ annihilation into neutrino (dav/d£i) or 
antincutrino (dav / '$,£2) of energy £1,2 measured in the 
plasma rest frame: 



dov 
d£\,2 



Gr 



8tt £i-E 2 



[(CV + C A f H h2 + (C v - C A ) 2 H 2A + 2m e 2 (C v 2 + C A 2 )H 3 ] G 



(31) 



where O is ( £,=1,2, £\ - v energy, £2 - v energy ) 



cated: 



and: 



O 



£± 




£t < £- 

£- < £i < £+ 

£i > £+ 



-(E 1 +E 2 )±-\p 1 + p 2 \ 



(32) 



(33) 



Hk = J2 h i\Pi+P2\ 



(34) 



With the index i = 1, 2 corresponding to neutrino and 
antincutrino, respectively, and 



Pi (P1+P2 cos 9) = 5, 
non-zero hi are: 



£i {E 1 +E 2 )-P 1 -P 2 -m e 



These values (£_,£+) restrict the neutrino energy £j 
to kinematically allowed (in a single e + e~ annihilation 
event) values. 

The neutrino cross-section is different from the an- 
tineutrino one, and the formula for neutrinos can be ob- 
tained from expression for antineutrinos by exchange of 
Hi and H 2 in (|3ip. as it is indicated by subscripts. 

Full expressions for Hi , H 2 and H3 in (|31|l are compli- 



hi 
hi 
hi 
h\ 
1'2 



2(A-£ l E 2 ) 2 

-4<5(A - £iE 2 )A - pfA 2 + £ 2 Pl 2 p 2 2 sin 2 1 

36 A 2 

2E\ 2 £ 2 

-A£iEx8A - pi 2 A 2 + £ 2 pi 2 p 2 2 sin 2 9 
3 5 A 2 

A-£ i (E 1 +E 2 ) 
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where Pi = |pj| denotes the length of 3-momentum and 
9 = Z(p l7 p 2 ) is the angle between incident electron and 
positron. 

The well-known formula for Dicus cross-section H17|) 
is reproduced by computing either of the following inte- 
grals: 



<TV D 



dav 



d£i, 



(36) 



The neutrino and antineutrino spectrum may be com- 
puted from the following formula: 

A(£i,a) =M J rf 3 P id 3 p 2 ^/i/a (37) 

where the normalization constant is related to the reac- 
tion rate, 



Af = R-\ 



(38) 



and can be computed from the formula 1)25(1 . Reaction 
rate can be computed from the cross-section (1311) as well, 
but due to lj3l)|l these expressions are equal. This is also 
obvious from physical point of view: the number of reac- 
tions is equal to number of neutrinos (or antineutrinos) 
emitted. 

The integral (|3T|) above, because of the presence of 
the unit step function in the cross-section (|31|l . must be 
evaluated numerically 4 . 

The most reliable method for multidimensional inte- 
gration is the Cuhre algorithm, but it is also possible 
to calculate integrals (|37|l using Monte Carlo algorithms. 
In actual calculations we used the Cuba library |24| . To 
get the complete neutrino spectrum, one has to compute 
numerically three-dimensional integral at every point. 

Sample results produced by our code p3| are given in 
Fig. |21 and in Fig. 0] Noteworthy, Monte Carlo integra- 
tion apparently fails to compute spectrum for neutrino 
energy £ < 0.1 MeV. However, taking into account pre- 
dicted purpose of calculations, these errors are insignifi- 
cant. Typically, spectrum will be used to estimate signal 
in neutrino detectors. Energy threshold for neutrino de- 
tection is usually much higher than 0.1 MeV, and this 
part of the spectrum is not detected at all. Therefore we 
may take advantage of the Monte Carlo integration per- 
formance, and compute spectrum much more than ten 
times faster. For theoretical considerations, we however 
recommend use of much reliable Cuhre algorithm. 

In Fig. 0] we can see clearly differences between v e , 
Vei ^/j,t and Vfir neutrino spectra. Particularly, electron 
anti-neutrino spectrum, main goal of e. g. GADZOOKS! 
detector has lower mean energy and tail than other 
neutrinos. 



Quality of the results computed using PSNS [23| may 
be judged using data from Table [I] where obtained nu- 
merically spectrum was integrated to get average neu- 
trino energy (columns 4-7). Effectively this means that 
4-dimcnsional integrals have been computed. Mean ener- 
gies were computed again using results from Sect. llTTDl 
as combination of Fermi-Dirac integrals l(21|l with accu- 
racy of at least 10 -7 (columns 8-11, only 4 digits are 
shown). No significant discrepancies were found. More- 
over, for fi <C kT or /i > kT, asymptotic expansions l|2"%|l 
and pnjl . respectively (columns 2-3), may be used as a 
very good estimate. 



D. Neutrino energy moments 

Fortunately, we were lucky to express moments of the 
neutrino spectrum by the Fermi-Dirac integrals l(21|l . 
With the moments given, one is able to approximate the 
spectrum with the aid of an appropriate analytical for- 
mula. 

Neutrino and antineutrino energy moments are com- 
puted as integrals 



£+ 

= J d£ h2 

£- 



d£ 



1,2 



(39) 



Unexpectedly, integration over the neutrino energy d£\ 
and the angle between p! and p 2 in l|39|l can be done 
analytically for any integer value 5 of n. However, we are 
unable to find general expression valid for any n similar 
to Eq. (HSl- 

For n = 0, due to (gBJ), Eq. (jSHJl reduces to Eq. 
with n = m = 0. Physically, it expresses the fact, that 
the reaction rate is equal to the number of emitted neu- 
trinos or antineutrinos per unit volume and per unit time. 

For n = 1, Eq. I|39|) is equal to the neutrino (antineu- 
trino) emissivity. For convenience, we express it in the 
form: 



Jl>* = Qx, 2 = | T AQ, 
where Q is the total emissivity (|24|l and AQ is: 



Q 



(40) 



Gp 
60TT 



4 G i~ G l 



+ r 



4 f G G_ 1/2 - G G 1/2 1 - ( G x G_ l/2 - G 1 G j 
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(g G\ 



/2 



G G 1 ^ 2 



1/2 



(41) 



From (|41|l we can see explicitly when neutrino and an- 
tineutrino emissivities (as well as the spectra) are not 
identical. This happens for 



4 Our code PSNS j^l will be available at 
http : //th-www . if . uj . edu . pl/psns 



5 At least up to n = 4 where we stopped calculations. 
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TABLE I: Comparison of the average neutrino energy {£) computed from numerically calculated (using 1371 implemented in 
psT ]) spectrum and formulae Chemical potential was chosen to be equal to 0.85 + m e MeV. Columns 2 and 3 show 

results obtained from asymptotic expansions I30H and 12811 . Up to the first order no differences between electron and mu/tau 
flavor exist. Columns 4-7 show mean neutrino energy as computed using spectrum from our PSNS |23| code employing CUHRE 
deterministic multidimensional integration algorithm from Cuba f24ll library. Columns 8-11 show results computed from our 
formulae 127H (mean neutrino-antineutrino pair energy) corrected using 1411 (mean neutrino energy is, of course Q/R where 
reaction rate is given by 1251 ) to get separately neutrino and antineutrino energies. Results are in perfect agreement (compare 
columns 4-7 and 8-11) even up to 4 digits. As Fermi-Dirac integrals 12 U (used in calculations presented in columns 8-11) are 
easily computed with accuracy of at least several digits |26|. we may conclude, that our code |23j| performs very well, and 
produce reliable results with accuracy possibly much better than 5% as estimated by Cuba library algorithms. Our analytical 
expansions 1301 128H are very useful (cf. columns 2-3 versus column 12 or 13) in the appropriate (degenerate or non-degenerate 
relativistic) conditions, as long as one do not need to know about differences between neutrino flavor, i.e neutrino-antineutrino 
asymmetry and e-[i,T asymmetry. 



kT [MeV] 
(/x = 1.361) 


lfx + 2kT 4W6kT 
+m e /2 


(£u.) <&.) <£„„. T > <&„, T ) 


(£„.> (£ De ) (£„ MiT ) (%, T ) 




0.01 


0.812 0.04 


0.703 0.913 0.753 0.911 


0.7037 0.9171 0.7532 0.9168 


0.810 0.835 


0.05 


0.900 0.205 


0.759 0.955 0.822 0.946 


0.7584 0.9544 0.8211 0.9469 


0.856 0.884 


0.1 


1.000 0.411 


0.852 1.029 0.930 1.023 


0.8510 1.028 0.9291 1.024 


0.940 0.977 


0.5 


1.800 2.053 


2.149 2.267 2.2495 2.2848 


2.1494 2.2672 2.2496 2.2854 


2.208 2.268 


1.0 


2.800 4.106 


4.1260 4.2329 4.2043 4.2333 


4.1277 4.2349 4.2050 4.2344 


4.181 4.220 


2.0 


4.800 8.212 


8.1912 8.2952 8.2507 8.2781 


8.1968 8.3004 8.2557 8.2833 


8.249 8.270 


5.0 


10.80 20.53 


20.476 20.578 20.525 20.551 


20.495 20.597 20.540 20.568 


20.546 20.554 


10.0 


10.80 41.06 


40.97 41.08 41.02 41.05 


41.019 41.121 41.053 41.081 


41.07 41.07 


100.0 


100.8 410.6 


410.1 410.2 410.1 410.1 


410.56 410.66 410.60 410.62 


410.61 410.61 



1. Parity violating interaction i.e. both Cy and Ca 
in © must be non-zero 

2. 7^ G~ i.e. distributions of electrons and 
positrons must be different. As the only difference 
between and G~ comes from chemical potential 
(cf. Eq. lj2*T|lL fi e must not be negligible. 

However, relative difference between neutrino and an- 
tineutrino emissivity usually is very small. For example, 
in the non-degenerate and relativistic case the total emis- 
sivity is: 

Q = 7 ~^- G F 2 (C v 2 + Ca 2 ) (kT) 9 (42) 

where 7£(5)/(127r) = 0.1925, while the difference between 
v and v luminosity is only: 

ja<j _(|£_«c^m) Gr , CvCAItm > 

(43) 

with the numerical coefficient in the parentheses equal to 
0.0043. 

We also present the second moment of the neutrino 
spectrum. Again, the formula is split into the average 
value {J 2) and a small additive term AJ 2 : 




0.01 0.1 1 10 

E [MeV] 

FIG. 3: Spectrum of the /i and r antineutrinos, computed 
from I37H with the use of the Cuhre algorithm implemented 
in Cuba library |2~| used in our PSNS code [23jl. Guaranteed 
relative accuracy is everywhere better than 3%. This figure 
explains why Cuhre deterministic algorithm (solid line), in 
spite of slow convergence, is recommended to compute neu- 
trino spectrum. Failure of Monte Carlo algorithms (dashed 
and dotted line) is apparent. However, these failures do not 
influence energy moments (mean energy and dispersion of the 
spectrum) significantly leading to errors smaller than those 
Monte Carlo algorithm produces itself. Detailed knowledge 
of the spectrum below 0.1 MeV seems also unimportant from 
experimental point of view. Temperature and chemical po- 
tential values as in Fig. |5] 



J 2 M = (J 2 ) ±AJ 2 . 



(44) 
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fitting formula / with n parameters 

A(£)~/(£,£l, &,...£„) 



(47) 



E [MeV] 

FIG. 4: On semilog plot differences between neutrino spectra 
of electron neutrinos (dashed), electron antineutrinos (solid), 
fi (r) neutrinos (dotted) and fi (t) antineutrinos (dash- 
dotted) are clearly visible. All four spectra are normalized 
to 1. T and fi the same as in Fig. [5] 



The average second moment is equal to: 



we may require that our fitting formula f(£) has exactly 
the same moments as the original spectrum A (£), given 
by e.g. E3: 



J? = J £?f(£i,Si,Z2, ■ ■ ■£„) d£i = J £?\(£i) d£i (48) 



Eqs. I|48[) form a set of algebraic equations with un- 
known parameters £i, £2, ■ ■ ■ £«■ 

Actually, only high energy tail of the neutrino spec- 
trum can be used to detect neutrinos 6 . To find the high 
energy behavior of neutrinos produced in the process it is 
not enough to know average v energy. We must compute 
at least the second moment. We face immediately the 
problem of appropriate explicit form for 1471) . For con- 
venience we want to use as simple formula as possible. 
Very successful approximation of the neutrino spectrum 
produced in various situations is given by the following 
formula W\: 



Gt 



777^ |(CV 2 +CU 2 ) 28 (g 3/2 G+ /2 + G+ /2 G 1/2 
30 G 1 Gf - 7 (g 3/2 G+ 1/2 + G+ /2 G_ 1/2 



£' 



(p 1/2^ -1/2 + G\i 2 G_ x i 2 



70 



( G 1/2 G -1 



-1/2 1 ^1/2^-1/2 j 



{s)a+ i ■«? 
Moments for this formula are: 

rxr-i 

(£) \"r(a + l + n) 



a + 1 

'HQ 1 

(49) 



) +60 (GrG+ + G+G -)] J £ n <t>{£) d£ = M (J&A 



+ 30 {Cv z - C A Z )G+G + 16 (3G y " - 2C A 2 )G~ /2 G+ /2 



(34G y 2 - 6G A 2 )GI 1/2 G+ 1 



1/2 



(45) 



r(a + l) 
u 

Particularly, the normalization constant A/" is: 

(a + l) a+1 



(50) 



A/": 



To obtain moments for neutrino or antineutrino spec- 
trum we have to add to (|45|l the term 1LAJ2: 



T(a + 1) ' 



(51) 



AJ, = 



G^ 2 
FStT 



CvCa 



4 I G 3 / 2 G x j 2 C 3/2 G 1 ^ 2 



the first moment is, of course, equal to the average v 
energy (£): 



(~<— s~t+ _ /~<+ 

^3/2^-1/2 ur 3/2 lT -l/2 



6 (G-^ Gq^ G-J^Gq 



^1/2^-1/2 ^1/2^-1/2 



(46) 



Higher neutrino energy moments could be computed 
also, but very good approximation for the whole spec- 
trum can be obtained from the two first moments. 



E. Fitting formula 



£ 4>{£) d£ = (£) (52) 
and dispersion is given by the following expression: 



(£ - (£)?4>(£) d£ 

\ J Va + 1 

N 



(53) 



Let's assume that the n-th moment J™ 2 of neutrino 
spectrum A(£) is known, for example from simula- 
tions, multidimensional integrations or numerical com- 
putations. Then, if we want to find parameters of the 



Analytical description of the high-energy tail of the spectrum 
could be very convenient, as the most of the existing big neutrino 
detectors operate only for £ > 4 MeV. Authors, however, have 
failed to find analytical formulae for the tail of the spectrum. 
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Dispersion ag is a rough measure how far from the mean 
energy the spectrum extends. 

If the spectrum is known, we can compute {£) and 
as which can be also computed from 1)2 7|l and (|45|l . One 
may also consider fitting the formula (|49|l to the spectrum 
in the least-squares sense if needed, but computing the 
moments is more straightforward. 

Comparison of the spectrum and the fit (|49|l is pre- 
sented in Fig. [SJ 

Given the spectrum from Monte Carlo simulations or 
computed using formula (|37[) . it is easy to find parameters 
a and (£) of the fitting formula l)49|l for any value of the 
chemical potential and temperature of the electron gas. 



IV. SUMMARY 

We have thoroughly analyzed details of the neutrino 
pair-annihilation process in the electron-positron plasma. 
Plasma is assumed to be in thermal equilibrium defined 
by the temperature T and the chemical potential [x. 
Given previously known results, based on the article of 
Dicus Q, we are able to compute combined neutrino- 
antincutrino emissivity Q and the mean energy (£) . The 
latter however, was not presented in the explicit form, 
so we derived the appropriate formula (|27|l . Some use- 
ful analytical expressions (|28II30|) in various regimes are 
also presented for the neutrino energy. Further progress, 
namely derivation of separate neutrino and antineutrino 
emissivities and spectra is impossible with the use of Di- 
cus cross-section (|17f) . Therefore we have derived cross- 
section (|31ll for pair annihilation into neutrino (antineu- 
trino) of given energy in the plasma rest-frame. Next 
step, derivation of the spectrum becomes quite simple. 
This spectrum is in full agreement (cf. Fig. |2J) with 
results of the Monte Carlo simulation. Unfortunately, 
the value at any single point of the spectrum requires 
evaluation of the resulting (effectively three-dimensional) 
integral H37[) . which can only be evaluated numerically 
due to the presence of the unit step function 1321) in 
the cross-section (|31|l . Therefore, we provide formulae 
(|24I41I45I46|) for neutrino spectrum moments as a com- 
bination of the Fermi-Dirac integrals. Using just two first 
moments, we are able to calculate parameters for the an- 
alytical approximation of the spectrum l|49[) . In some 
cases, this procedure gives particularly simple expres- 
sions. For example, in relativistic and non-degenerate 
regime (kT > 2m e , kT > /i) the spectrum is given by: 




FIG. 5: Numerically computed points of the neutrino spec- 
trum (symbols) and the fit (lines). Curves represent spectra 
for kT = 1, 5, 10 MeV. Even in in the worst case considered, 
kT =1 MeV (o), fit i^4l is particularly good. 



A = 



(a + l) a+1 /2700C(5) 

T(a + 1) v 771-4 



(57) 



Numerically: a = 3.180657028, A = 0.1425776426, 
a = 1.018192299. Previously, the spectrum given by the 
analytical formula (|54|l above, must have been computed 
by the means of the Monte Carlo simulation! Compari- 
son of the spectrum computed numerically and given by 
eq. (|54|l is presented in Fig. |SJ 

In more difficult cases, e.g. required for neutrinos from 
the pre-supernova core, where kT ~ m e ~ fi, one must 
use general expressions for neutrino energy moments, 
with wealth of analytical [T^ . |2S| and numerical [2{J [2|J 
methods available for calculating Fermi-Dirac integrals. 
In the degenerate case differences between neutrino and 
antineutrino spectrum (and between v e and v^ ^ as well, 
cf. Fig. 0} must be taken into account. However, in 
typical situations they may be considered as a small per- 
turbation to the average value. 

If for some reasons exact results are required, the spec- 
trum may be computed from i|37|) point-by-point as well. 



A 



kT \kT 



with parameters: 



£ 



exp (—a£/kT) 



56700 C(5) 7T 4 



217 7T 10 - 13668750 C(5) 2 

217tt 10 - 35538750 C(5) 2 
217tt 10 - 13668750 C(5) 2 



(54) 

(55) 
(56) 
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